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^"*^ ' We give a new and detailed description of the structure of cut loci, with direct applications to the singular 

sets of some Hamilton- Jacobi equations. These sets may be non-triangulable, but a local description at all 
points except for a set of Hausdorff dimension n — 2 is well known. We go further in this direction by giving a 
clasification of all points up to a set of Hausdorff dimension n — 3. 
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Abstract 



1 Introduction 



In this paper we improve the current knowledge about the sets known as the cut locus in differential geometry and 
! the singular set of solutions to static Hamilton- Jacobi equations: 

H(x,du(x)) = l xeM (1.1) 

u(x) = g(x) x e dM (1.2) 

for H smooth and convex in the second argument and g satisfying a standard compatibility condition (see !3.ip 
' The solution to the equations above is given by the Lax-Oleinik formula: 

(N ' 

(N ■ «(p)= {d(p,q) + g(q)} (1.3) 

• ■ q£OM 

where d is the distance function of a Finsler metric constructed in fl from the hamiltonian function H . Thus, when 
.9 = 0, the solution to the equations is the distance to the boundary, and then the singular set of the solution is the 
cut locus from the boundary (see |LNj ). an object of differential geometry. In section[3]we find a similar relationship 
when g^O. 

Our main result is a local description around any point of the cut locus except for a set of Hausdorff dimension 
n-3 (see Theorem |2~2|) . 

This structure result was originally motivated by its use in the paper [AG] , This application motivated some 
important decisions. For example, all the proofs apply to the more general balanced split locus. We show in this 
paper that cut loci (hence singular sets of solutions to HJ equations) are balanced split loci. In general, there are 
many balanced split loci besides the cut locus. In the paper |AG| . and using the results in this paper, we study and 
classify all possible balanced split loci. 

We believe that our description of the cut locus could also be useful in other contexts. For instance, the study 
of brownian motion on manifolds is often studied on the complement of the cut locus from a point, and then the 
results have to be adapted to take care of the situation when the brownian motion hits the cut locus. As brownian 
motion almost never hits a set with null H n ~ 2 measure, we think our result might be useful in that field. 

The paper is divided in six sections besides this introduction and an appendix. For the convenience of the reader 
we have included separate statements of our results in section [2] together with examples showing that some of them 
are sharp, a compendium of previous results in the literature and suggestions for future work. In section [3] we 
enlarge the class of Hamilton-Jacobi problems for which our results apply: this allows to expand the applicability of 
a result by Li and Nirenberg (cf. [LN]). Section [4] contains all the necessary definitions that we use along the paper; 
although some of them have already appeared elsewhere, we have considered useful to collect them here in order to 
save the reader some effort. More important, this section contains also the key notions of split locus and balanced 
split locus, that play a key role in the rest of the paper. In section [5] we show that the cut locus of a submanifold 
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in a Finsler metric is a balanced set. This is an extension of the corresponding Riemannian claim in [IT2] , and it 
is necessary in order to apply our results in situations requiring the extra Finsler generality, as for instance in the 
already mentioned Hamilton- Jacobi problems. Section [6] proves our results concerning focal vectors in a balanced 
split locus (in the context of a cut locus, focal minimizing geodesies), and section [7] contains the results about the 
structure of balanced split loci up to codimension 3. An Appendix contains some important facts about Finsler 
exponential maps. 

Acknowlegdements The first author came upon this problem after working with Yanyan Li, who gave many 
insights. The authors benefited from conversations with Luc Nguyen and Juan Carlos Alvarez Paiva. Both authors 
were partially supported during the preparation of this work by grants MTM2007-61982 and MTM2008-02686 of 
the MEC and the MCINN respectively. 

2 Statements of results 

2.1 Setting 

From now on, we will work in the following setting: 

• A C°° Finsler manifold M with compact boundary dM. The space M U dM need not be compact. 

• The geodesic vector field r in TM. 

• A smooth map M : dM — * TM that is a section of the projection map it : TM — » M of the tangent to M, 
and such that T(x) points to the inside of M for every x £ dM. 

Let $ be the flow of r, and D($) its domain. We introduce the set V: 

V = {#(*, r(as)), t > 0, x £ dM, (t, T(x)) e £>(*)} (2.1) 

The interior of V is locally invariant under $t (equivalently, r is tangent to V). We set F to be the map ir\v ■ V — > M. 

We say a point x e V is a focal point iff d x F is a singular map, and call dimker(d x F) the order of x. Finally, 
let S be a balanced split locus for this setting. 

Remark. Our results covers both the cut locus from a point and the cut locus from a hypersurface. However, 
let us recall that, when the interest is in the cut locus, we only need to consider the exponential map from an 
hypersurface. The cut locus of a point is also the cut locus of a small sphere centered at the point. In this way, 
our focal points with respect to the sphere are the conjugate points with respect to the point. The cut locus of a 
smooth submanifold is also the cut locus of an e-neighborhood of the submanifold. 

Observe also that some authors use the term conjugate instead of focal, even when studying the distance function 
from a hypersurface (see for instance |LN] ). 

2.2 Results 

We will show that a cut locus is a balanced split locus (see section |4] for the definition of this term and section [5] 
for the proof), so the reader may simply think that the following results apply to the cut locus. In this situation, 
the set Rp with p 6 M consists of the vectors tangent to the minimizing geodesies from p to dM. Nonetheless, the 
notation for the general case is explained in definition 14.41 

Our main result asserts that we can avoid focal points of order 2 and above if we neglect a set of Hausdorff 
dimension n — 3.: 

Theorem 2.1 (Focal points of order 2). There is a set N C S of Hausdorff dimension at most n — 3 such that for 
any p € S \N and x £ V such that F(x) = p and d x F{r x ) G R p : 

dim(kerd 2; F) < 1 

Combining this new result with previous ones in the literature, we are able to provide the following description 
of a cut locus. All the extra results required for the proof of these result will be proved in this paper, for the 
convenience of the reader, and also because some of them had to be slightly generalized to serve our purposes. 
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Theorem 2.2 (The cut locus up to 7i-codimension 3). Let S be either the cut locus of a point or submanifold in a 
Finsler manifold or the closure of the singular locus of a solution of \l.l\ and \1.2\ Then S consists of the following 
types of points : 

• Cleave points: Points at which R p consists of two non-focal vectors. The set of cleave points is a smooth 
hypersurface; 

• Edge points: Points at which R p consists of exactly one vector of order 1. This is a set of Hausdorff 
dimension at most n — 2; 

• Degenerate cleave points: Points at which R p consists of two vectors, such that one of them is conjugate 
of order 1, and the other may be non-conjugate or conjugate of order 1. This is a set of Hausdorff dimension 
at most n — 2; 

• Crossing points: Points at which R p consists of non-focal and focal vectors of order 1, and R* is contained 
in an affine subspace of dimension 2. This is a rectifiable set of dimension at most n — 2; 

• Remainder: A set of Hausdorff dimension at most n — 3; 

Finally, in regard to singular sets of viscosity solutions to HJ equations, we prove the following extension Theorem 
1.1 of |LN| . In this result dM may not be compact. 

Theorem 2.3. Let S be the singular set of a solution to the Hamilton- Jacobi system 

H(x,du(x))=l x e M 

u(x) — g(x) x £ dM 

where g : dM — > R is a positive smooth function such that \g(y) — g(z)\ < kd(y,z) for some k < 1. If fi is the 
function whose value at y G dM is the distance to S along the unique characteristic departing from y, then 

1. fi is Lipschitz. 

2. If in addition dM is compact, then the (n — 1)- dimensional Hausdorff measure of S n K is finite for any 
compact K . 

3. S is a Finsler cut locus from the boundary of some Finsler manifold. 
2.3 Examples 

We provide examples of Riemannian manifolds and exponential maps which illustrate our results. 

First, consider a solid ellipsoid with two equal semiaxis and a third larger one. This is a 3D manifold with 
boundary, and the geodesies starting at the two points that lie further away from the center have a first focal 
of order 2 while remaining minimizing up to that point. This example shows that our bound on the Hausdorff 
dimension of the points in the cut locus with a minimizing geodesic of order 2 cannot be improved. 

Second, consider the surface of an ellipsoid with three different semiaxis (or any generic surface as in [B] , with 
metric close to the standard sphere) and an arbitrary point in it. It is known that in the tangent space the set of 
first focal points is a closed curve C bounding the origin, and at most of these points the kernel of the exponential 
map is transversal to the curve C. More explicitely, the set C* of points of C where it is not transversal is finite. 
Consider then the product M of two such ellipsoids. The exponential map onto M has a focal point of order 2 at 
any point in (C\C*) x (C\C*), and the kernel of the exponential map is transversal to the tangent to C x C. 
Thus the image of the set of focal points of order 2 is a smooth manifold of codimension 2. 

This example shows that the statement of theorem 12.11 cannot be simplified to say only that the image of the 
focal points of order 2 has Hausdorff dimension n — 3. 

Finally, recall the construction in |GS] . where the authors build a riemannian surface whose cut locus is not 
triangulable. Their example shows that the set of points with a focal minimizing geodesic can have infinite H n ~ 2 
measure. A similar construction replacing the circle in their construction with a 3d ball shows that the set of points 
with a minimizing geodesic focal of order 2 can have infinite H n ~ 3 measure. 
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2.4 Relation to previous results in the literature 



Our structure theorem generalizes a standard result that has been proven several times by mathematicians from 
different fields (see for example |BL], (H], [MM] and [IT]): 

A cut locus in a Riemannian manifold is the union of a smooth (n — l)-dimensional manifold C and a set 
of zero (n — l)-dimensional Hausdorff measure (actually, a set of Hausdorff dimension at most n — 2). 
The set C consists of cleave points, which are joined to the origin or initial submanifold by exactly two 
minimizing geodesies, both of which are non-focal. 

We observe that this theorem follows from our theorem 12,21 since the union of edge, degenerate cleave, and 
crossing points is a set of Hausdorff dimension at most n — 2. Our main contribution is to show that, up to 
codimension 3, these latter ones are the only new type of points that can appear. 

The statement on cleave points quoted above follows from lemmas 17.21 17.31 an d 16.31 only. Theorem 12.11 is not 
necessary if a description is needed only up to codimension 2. The proof of the three lemmas is simple and has 
many features in common with earlier results on the cut locus. However, we have decided to include a proof of 
them that applies to balanced split loci, because not every balanced split locus coincides with the cut locus (see 
[AG]), and the extra generality is necessary for forthcoming work. 

In a previous paper, A. C. Mennucci studied the singular set of solutions to the HJ equations with only C k 
regularity. Under this hypothesis, the set S\C may have Hausdorff dimension strictly between n — 1 and n — 2 
(see [Mj). We work only in a C°° setting, and under this stronger condition, the set S \ C has always Haussdorf 
dimension at most n — 2. 

Our result [22] uses the theory of singularities of semi-concave functions that can be found for example in [A AC] , 
Though their result can be applied to a Finsler manifold, we had to give a new proof that applies to balanced sets 
instead of just the cut locus. 

Finally, the very definition of balanced split locus is inspired in lemma 2.1 of [IT2]. Slight changes were required 
to adapt the property to Finsler manifolds, and the proof of the lemma itself. 

2.5 Further questions 

Theorem 12.11 and the classical result quoted earlier suggest the following conjecture: although the image of the 
focal points of order k in an exponential map can have Hausdorff dimension n — k, the set of points in M with a 
minimizing geodesic of order k only has Hausdorff dimension n — k — 1. 

The examples in the above section can be extended to focal points of greater order without pain, showing that 
this conjecture cannot be improved. 

In this paper all the structure results about cut loci follow from the split and balanced properties of a cut locus. 
We will address the question of how many balanced split sets are there in a future paper. We believe this approach 
is an interesting way to look at viscosity solutions and their relation with classical solutions by characteristics. 

Finally, we would like to mention that similar hypothesis and similar structure results hold in other settings. 
It would be interesting to study the structure of the singular locus of the solutions to other Hamilton-Jacobi 
equations, when the Hamiltonian depends not only on x and du, but also in t and u itself, for the Dirichlet and 
Cauchy problems, or maybe without the convexity hypothesis on H. 

3 Singular locus of Hamilton-Jacobi equations 

In this section we study the relationship between Hamilton-Jacobi equations and Finsler geometry. The reader can 
find more details in |LN] and [L]. 

Let M be an open set (or manifold) M with possibly non-compact boundary. We are interested on solutions to 
the system 



where H : T*M — > M is a smooth function that is 1-homogeneous and subadditive for linear combinations of 
covectors lying over the same point p, and g : dM — > R is a smooth function that satisfies the following compatibility 
condition: 



H(x, du(x)) = 1 
u(x) = g(x) 



x £ M 



x 6 dM 



\g(y)-g(z)\<kd(y,z) Vy,zedM 



(3.1) 
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for some k < 1, 

As is well known, the unique viscosity solution is given by the Lax-Oleinik formula: 

u(p)= inf ■ {d(p,<z) + fl(g)} (3.2) 

where d is the distance induced by the Finsler metric that is the pointwise dual of the metric in T*M given by H: 

(p p (v) = sup{(w,a) p : a € IJM, H(p,a) = l} (3.3) 

A local classical solution can be computed near dM following characteristic curves, which are geodesies of the 
metric (p starting from a point in dM with initial speed given by a vector field on dM that we call the characteristic 
vector field. The viscosity solution can be thought of as a way to extend the classical solution to the whole M . 

When .g = 0, the solution (|1,3|1 is the distance to the boundary. It can be found in [LN] . among others, that the 
closure of the singular set of this function is the cut locus, given for example by: 

g _ f x ]\fl. there are at least two minimizing geodesies from dM to x 1 ^ 
\ ' or the unique minimizing geodesic is focal J 

Hamilton-Jacobi equations fit our setting if we let the vector field r be the geodesic vector field, and T be the 
vector field at dM that is tangent to the departing characteristics. The map F : V — > M is the map sending 
(x,t) S dM x 1 to lv(x){t)i f° r the geodesic 7 with initial speed v(x), where v : dM — > SM is the characteristic 
vector field, and V C dM x R is the domain of definition of F. The characteristic vector at x is the inner pointing 
normal if g = (see the appendix for the definition of normal under Finsler conditions). 

Our intention in this section is to adapt this result to the case g > 0. If dM is compact, a global constant 
can be added to an arbitrary g so that this is satisfied and S is unchanged. We still require that g satisfies the 
compatibility condition 13.11 Under these conditions, our strategy will be to show that the Finsler manifold (M,(p) 
can be embedded in a new manifold with boundary (N, (p) such that u is the restriction of the unique solution u to 
the problem 

H(x,du(x)) = l xEN 
u(x) =0 x edN 

thus reducing to the original problem (H and (p are dual to one another as in I3.3fl . This allows us to characterize 
the singular set of (|1.3p as a cut locus, as well as draw conclusions similar to those in [LN] . 



Definition 3.1. The indicatrix of a Finsler metric tp at the point p is the set 

I P = {veT p M : <p(p,v) = l} 
Lemma 3.2. Let tpa and (pi be two Finsler metrics in an open set U, and let X be a vector field in U such that: 

• The integral curves of X are geodesies for ipo . 

• <Po(p,X P ) = <Pi(p,X P ) = 1 

• At every p £ U , the tangent hyperplanes to the indicatrices of ipo and (pi in T p U coincide. 
Then the integral curves of X are also geodesies for (p\ 

Proof. Let p be a point in U. Take bundle coordinates of T P U around p such that X is one of the vertical coordinate 
vectors. An integral curve a of X satifies: 

(cp Q ) p {a(t),a'(t)) = (<p 1 ) p (a(t),a'(t)) = 

because of the second hypothesis. The third hypothesis imply: 

{(po) v (a(t),a'{t)) - (vi) v (a(t),a'(t)) 

So inspection of the geodesic equation: 

(p p (a(t),a'(t)) = j t ((p v (a{t),a'(t))) (3.5) 
shows that a is a geodesic for (p\, □ 
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Corollary 3.3. Let ip be a Finsler metric and X a vector field whose integral curves are geodesies. Then there is 
a Riemannian metric for which those curves are also geodesies. 

Proof. The Riemannian metric gij(p) = d ® v . (pip, X) is related to <p as in the preceeding lemma. □ 

Lemma 3.4. Let X be a non-zero geodesic vector field in a Finsler manifold and u> its dual differential one-form. 
Then the integral curves of X are geodesies if and only if the Lie derivative of u> in the direction of X vanishes. 

Proof. Use lemma l3"T3l to replace the Finsler metric with a Riemann metric for which w is the standard dual one-form 
of X in Riemannian geometry. Now the lemma is standard. □ 

Proposition 3.5. Let M be an open manifold with smooth boundary and a Finsler metric (p. Let X be a smooth 
transversal vector field in dM pointing inwards (resp. outwards). Then M is contained in a larger open manifold 
admitting a smooth extension (p of ip to this open set such that the geodesies starting at points p € dM with initial 
vectors X p can be continued indefinitely backward (resp. forward) without intersecting each other. 

Proof. We will only complete the proof for a compact open set M and inward pointing vector X, as the other cases 
require only minor modifications. 

We start with a naive extension ip' of tp to a larger open set M 2 D M. The geodesies with initial speed X can 
be continued backwards to M 2 , and there is a small e for which the geodesies starting at dM do not intersect each 
other for negative values of time before the parameter reaches —e. 

Define 



P : dM X (-e, 0] -»■ M 2 , P(q,t) := a q {t) 

where a q : (—6,0} — > M 2 is the geodesic of ip 1 starting at the point q € dM with initial vector X q . When 
p G U e := Image(P) there is a unique value of t such that p = P(q,t) for some q £ dM. We will denote such t by 
d{p). Extend also the vector X to U e as X p = ci q (t) where p = P(q,t). 
Let c : (—£,0] — > [0, 1] be a smooth function such that 

• c is non-decreasing 

• c(t) = 1 for - e/3 < t 

• c(t) = for t < -2e/3 
and finally define 

X p = c(dip))X p 

in the set U e . 

Let wo be the dual one form of X with respect to (p for points in dM, and let u> be the one form in U e whose Lie 
derivative in the direction X is zero and which coincides with loq in dM. Then we take any metric ip" in U s (which 
can be chosen Riemannian) such that X has unit norm and the kernel of oj is tangent to the indicatrix at X. 

By lemma I3T41 the integral curves of X are geodesies for ip". Now let p be a smooth function in U e U M such 
that p\m — lj p\u E \u e/3 — and < p < 1, and define the metric: 

<P = P(P)<P(P> v) + (l- pip))(p"ip, v) 

This metric extends ip to the open set U s and makes the integral curves of X geodesies. As the integral curves 
of X do not intersect for small t, the integral curves of X reach infinite length before they approach dll e and the 
last part of the statement follows. □ 

Application of this proposition to M and the characteristic, inwards-pointing vector field v yields a new manifold 
N containing M, and a metric for N that extends (p (so we keep the same letter) such that the geodesies departing 
from dM which correspond to the characteristic curves continue indefinitely backwards without intersecting. 

This allows the definition of 

P:dMx (-co, 0] -> N, P(q, t) := a q (t) 
where a are the geodesies with initial condition X, continued backwards. Finally, define u : U — » R by: 

fi(s) = * = ^ N \ M (3.6) 

I u(x) x e M 



6 



We notice that both definitions agree in an inner neighborhood of dM, so the function u is a smooth extension of 
u to N. 

Theorem 3.6. Let A = u~ 1 (0). Then the following identity holds in {u > 0} : 

u(x) = d(x, A) (3.7) 
Proof. Let g t be the flux associated to the characteristic vector field X. By definition of u, we see that: 

g*u(x) = u(x) + t 

at least for (x,t) in an open set O containing N\M X (— oo, 0]. We deduce that gt, restricted to a small ball B, 
sends the intersection of a level set of u with the ball to another level set of u, whenever t is small enough so that 
gt(B) is contained in O. 

In particular, the tangent distribution to the level sets is transported to itself by the flow of X . On the other 
hand, the orthogonal distribution to X is also parallel, so if we show that they coincide near dM, we will learn that 
they coincide in O. 

Now recall that inside M, u coincides with u, which is also given by the Lax-Oleinik formula fOl Let y G dM 
and t > small. This formula yields the same value as the local solution by characteristics, and we learn that the 
point y is the closest point to gt(y) on the level curve {u — u(y)}. By appeal to lemma 2.3 in |LN| . or reduction to 
the Riemannian case as in l3.3[ we see that the level set {u = uq} is orthogonal to the vector X y . It follows that, in 
O: 

H(x,du(x)) = sup{du(x)(Y) : ip(Y) = 1 } = du(x)(X) = 1 

In order to show that u and d( ■ , A) agree in U, we use the uniqueness properties of viscosity solutions. Let N 
be the open set where u > 0. The distance function to A is characterized as the unique viscosity solution to: 

• u — in A 

• H(x,du(x)) = 1 in N 

Clearly u satisfies the first condition. It also satisfies the second for points in the set M because it coincides 
with u, and for points in N \ M because H(x, du{x)) — 1 there. 

□ 

Proof of Theorem \2.3l The first part follows immediately from Theorem 13.61 and Theorem 1.1 in |LN] , The second 
is an easy consequence of the first, while the last is contained in the results of this section. □ 

Remark. Regularity hypothesis can be softened. In order to apply the results in |LN| . it is enough that the 
geodesic flow, the characteristic vector field and g itself are C 2,1 , which implies that A is C 2,1 . Thus the result in 
true for less regular hamiltonians and open sets. 

4 Split locus and balanced split locus 

We now introduce some properties of a set necessary in the proofs of our results. We prove in section [5] that a cut 
loci in Finsler manifolds have all of them. 

Definition 4.1. For a pair of points p,q G M such that q belongs to a convex neighborhood ofp, we define, following 

v P (q)=7(0) (4.1) 
as the speed at of the unique unit speed minimizing geodesic 7 from p to q. 

Definition 4.2. The approximate tangent cone to a subset E at p is: 

T(E,p) = {r6:9 = \imv p { Pn ), 3{p n } C E, Pn ^p,r> 0} 
and the approximate tangent space Tan(E,p) to E at p is the vector space generated by T(E,p). 

We remark that the definition is independent of the Finsler metric, despite its apparent dependence on the 
vectors v p ( Pn ). 
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Definition 4.3. For a set S C M , let A(S) be the union of all integral segments of r with initial point in T whose 
projections in M do not meet S. We say that a set S C M splits M iff ir restricts to a bijection between A(S) and 
M\S. 

Whenever S splits M, we can define a vector field R p in M \ S to be dF x (r x ) for the unique x in V such that 
F{x) — p and there is an integral segment of r with initial point in T and end point in x that does not meet F~ 1 (S). 

Definition 4.4. For a point p £ S , we define the limit set R p as the set of vectors in T p M that are limits of 
sequences of the vectors R q defined above at points q £ M \ S. 




Figure 1: An arbitrary split locus and a balanced split locus 



Definition 4.5. A set S that splits M is a split locus iff 

S = {peS: $R P > 2} 

The role of this condition is to restrict S to its essential part. A set that merely splits M could be too big: 
actually M itself splits M. Finally, we introduce the following more restrictive condition. 

Definition 4.6 (Balanced split locus). We say a split locus S C M is balanced at p £ S iff for any sequence {p n } 
converging to p with v Pn (p) and X n € R Pn approaching v £ T X M and £ R p respectively, then 

Woo{v) = max{w(ti) : w is dual to some R £ R p } 

where Wqo is the dual of Xoq. We say S is balanced if it is balanced at every point. 



5 Balanced property of the Finsler cut locus 

In this section we show that the cut locus of a Finsler exponential map is a balanced set. The proof is the same as 
in lemma 2.1 in |IT2j . only adapted to Finsler manifolds, where angles are not defined. 

Proposition 5.1. The cut locus of a Finsler manifold M with boundary is a balanced split locus. Moreover, for p, 
p n , v and Xoo as in the definition of a balanced split locus, we have 

v d(dM,p)-d(dM,p n ) , , 

hm — = Wooiv) 

n^oo d(p,p n ) 
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Proof. The cut locus S splits M, as follows from the well-known property that if a geodesic 7 from dM to p = j(t) 
is minimizing, and s < t, then 7|[o, s ] ls the unique minimizing geodesic from dM to 7(s), and is non-focal. 

It is also a split locus, as follows from the characterization of the cut locus as the closure of the singular set of 
the function distance to the boundary (as found in |LN] for example) . The distance to the boundary is differentiable 
at a point if and only if there is a unique minimizing geodesic from the point to the boundary. 

Next we show that S is balanced. Take any Y £ R p , and let 7 be the minimizing geodesic segment joining dM 
to p with speed Y at p. Take any point q 6 7 that lies in a convex neighborhood of p and use the triangle inequality 
to get: 

d(dM,p) - d(dM,p n ) > d(q,p) - d(q,p n ) 
Then the first variation formula yields, for a constant C: 

d(q,p) - d(q,p n ) > w(v Pn (p))d{p n ,p) - Cd(p,p n ) 2 

and we get: 

lhjiM d(dM,p) '-d(3M, Pn)> 

rc->oo d{p,p n ) 

for any w that is dual to a vector in R p . 

Then consider Xoo, let 7 be the minimizing geodesic segment joining dM to p with speed X^ at p, and let j n 
be the minimizing geodesic segment joining dM to p n with speed X n at p n , Take points q n in j n that lie in a fix 
convex neighborhood of p. Again: 

d(dM,p) - d(dM,p n ) < d(q n ,p) - d(q n ,p n ) 
while the first variation formula yields, for a constant C: 

d(q n ,p) - d(q n ,Pn) < w(v Pn (p))d(p„,p) - Cd(p,p n ) 2 

and thus: 

r d(dM,p)-d{dM,p n ) 

hmsup — < Woo (X) 

n^oo d(p, p n ) 

This proves the claim that S is balanced. □ 



6 Focal points in a balanced split locus 

In this section we prove Theorem 12.11 Throughout this section, M, r, V and F are as in section 12.11 and S is a 
balanced split locus as defined in 14.61 

Definition 6.1. A singular point x G V of the map F is an A2 point if ker(dF x ) has dimension 1 and is transversal 
to the tangent to the set of focal vectors. 

Remark. Warner shows in [Wj that the set of focal points of order 1 is a smooth (open) hypersurface inside V, 
and that for adequate coordinate functions in V and M, the exponential has the following normal form around any 
A2 point, 

(xi,X2, ■ ■ -,x m ) — ► {x\,x 2 , ■ ■ -,x m ) (6.1) 
Proposition 6.2. For any p E M and X £ R p , the vector X is not of the form dF x (r) for any A2 point x. 

Proof. The proof is by contradiction: let p £ S be such that R p contains an A2 vector Z. There is a unique c eV 
such that F(c) = p and dF c (r c ) = Z. By the normal form l|6.ip . we see there is a neighborhood U of c such that no 
other point in U maps to p. Furthermore, in a neighborhood B of p the image of the focal vectors is a hypersurface 
H such that all points at one side (call it Bi) have two preimages of F\u, all points at the other side B2 of H have 
no preimages, and points at H have one preimage, whose corresponding vector is A2-focal. It follows that Z is 
isolated in R p . 

We notice there is a sequence of points p n — > p in B 2 with vectors Y n e R Pn such that Y n — ► Y ^ X. Thus R a 
does not reduce to Z. 
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The vector Z is tangent to if, so we can find a sequence of points p n £ B 2 approaching p such that 



lim v Pn (p) = Z 

n — >oo n 

We can find a subsequence p nk of the p n and vectors X k £ -Rp„ fc such that converges to some Xoo £ R p . By the 
above, X x is different from Z, but Z(X) < 1 = Z(Z) (where Z is the dual form to Z), so the balanced property is 
violated. □ 

The following is the analogous to theorem 12,11 for focal points of order 1. 

Proposition 6.3 (Focal points of order 1). There is a set N C S of Hausdorff dimension n — 2 such that for all 
p 6 S \N and x € V such that F{x) — p and d x F(r x ) € R p , the linear map d x F is non-singular. 



Proof. The proof is identical to the proof of lemma 2 in [IT] for a cut locus, but we include it here for completeness. 
First of all, at the set of focal vectors of order k > 2 we can apply directly the Morse-Sard-Federer theorem (see 
[F]) to show that the image of the set of focal cut vectors of order k > 2 has Hausdorff dimension at most n — 2. 

Let Q be the set of focal vectors of order 1 (recall it is a smooth hypersurface in V) . Let G be the set of focal 
vectors such that the kernel of dF is tangent to the focal locus. Apply the Morse-Sard-Federer theorem again to 
the map F\q to show that the image of G has Hausdorff dimension at most n — 2. Finally, the previous result takes 
cares of the A2 points. □ 



We now turn to the main result of this paper: we state and prove Theorem 16.41 which has 12.11 as a direct 
consequence. In order to study the map F more comfortably, we define the special coordinates in a neighborhood 
of a focal point z of order k. 



Special coordinates. Let B = {v\, . . . ,v n } be the basis of T Z V indicated in the second part of Proposition 18. 3| 
and B' F r z -, the corresponding basis at F(z) e M formed by vectors d z F(vi), d\F{v\ (jua ),■■■, d%F{v\$Vk+i), and 
d z F(vi), i > k + 1. 

Make a linear change of coordinates in a neighborhood of F(z) taking &'f< z \ to the canonical basis. The 
coordinate functions F % (x) — F l (z) of F for i ^ 2, . . . , k + 1 can be extended to a coordinate system near z with 
the help of k functions having V2, ■ ■ ■ , «fc+i as their respective gradients at z. In this coordinates F looks: 



F(x 1 ,x 2 , ■ ■ - ,x k +i, ■ ■ ■ ,x n ) = (x 1 ,F z (x), 



F 



fc+1 



(x),X k+2 , ■ ■ ■ ,X n ) 



(6.2) 



Theorem 6.4. Let M, V, F and r be as in section POl Let S be a balanced split locus H4.6\ ). The set of focal 
points of order 2 in V decomposes as the union of two subsets Q\ and Q\ such that: 

• No vector in Q\ maps under dF to a vector in any of the R a . 

• The image under F of Q\ has Hausdorff dimension at most n — 3. 

Proof. Let z be a focal point of order 2 and take special coordinates at U z near z. In the special coordinates near 
z, F is written: 

..,x n ) = (x!,F z (x),F z (x),Xi, ...,x n ) (6.3) 
, , x n ) in a neighborhood U z of z with F(0, . . . , 0) = (0, . . . , 0). 



F(x 1 ,x 2 ,x 3 ,Xi, 
for some functions F z and F z , and x = {x\, . 



The Jacobian of F is: 



JF 



* 








3X2 


dx 2 


d£l 


dFl 




dx3 


* 


* 


* 





A point x is of second order if and only if the 2x2 submatrix for the x 2 and X3 variables vanish: 



VdFi, \ 



dx 2 W 



2£l( x ) 211 t x ) 



= 



(6.4) 
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We write: 

F 2 (x) = x x x 2 + q(x 2 ,x 3 ) + T 2 (x) 

F s z {x) = x lX3 + r(x 2 ,x 3 ) + T 3 (x) 

where q(x 2 ,x 3 ) and r(x 2 , x 3 ) are the quadratic terms in x 2 and x 3 in a Taylor expansion, and T consists of terms 
of order > 3 in x 2 and x 3 , and terms of order > 2 with at least one X4, i ^ 2, 3. x The nature of the polynomials q 
and r in the special coordinates at z will determine whether z is in Q\ or in Q 2 . We have the following possibilities: 

1. either q or r is a sum of squares of homogeneous linear functions in x 2 and x 3 (possibly with a global minus 
sign) . 

2. both q and r are products of distinct linear functional (equivalently, they are difference of squares). Later 
on, we will split this class further into three types: 2a, 2b and 2c. 

3. one of q and r is zero, the other is not. 

4. both q and r are zero. 

We set Q\ to be the points of type 1 and 2c, and Q\ to be the points of type 2a, 3 and 4. Points of type 2b do 
not appear under the hypothesis of this theorem. 



Type 1. The proof is similar to Proposition 16.21 Assume z = (0, . . . , 0) is of type 1. If, say, q is a sum of squares, 
then in the subspace given by x\ — a and x± = ■ ■ ■ = x n = 0, x 2 will reach a minimum value that will be greater 
than —Ca 2 for some C > 0. We learn there is a sequence p k — (t k ,—(C + l)(t k ) 2 , 0, . . . , 0), for t k f 0, approaching 
(0, . . . , 0) with incoming speed (1, 0, . . . , 0) and staying in the interior of the complement of F(U) for k large enough. 
Pick up any vectors Vk € R p k converging to some Vb (passing to a subsequence if necessary) . Then Vq is different 
from (1, 0, . . . , 0) € Ro, and 

v^((i,...,o))<(C77o)((i,...,o)) = i 

violating the balanced condition. 



Type 2 and 3. We take special coordinates at a fixed x and assume q =/= 0. Before we start, we will change 
coordinates to simplify the expression of F further. Consider a linear change of coordinates near x that mix only 
the x 2 and x 3 coordinates. 



1 x' 2 < 


).A. 


ill) 


\ 4 , 







followed by the linear change of coordinates near p that mix only the y 2 and y 3 coordinates with the inverse of the 
matrix above: 

' 2/2 



2/2 



= A~ x ■ 



27.3 



Straightforward but tedious calculations show that there is a matrix A such that the map F has the following 
expression in the coordinates above: 

F(xi,x 2 ,x 3 ,X4, ...,x n ) = (xi,xix 2 + (xj - xl),xix 3 + r(x 2l x 3 ),Xi, x n ) + T 

In other words, we can assume q{x 2 ,x 3 ) — [x\ — x 3 ). 

Take xi, fixed and small, i > 3. At the origin, JF is a diagonal matrix with zeros in the positions (2,2) and 
(3, 3). We recall that z is focal of order 2 iff the submatrix l|6.4p vanishes. This submatrix is the sum of 



Xi + 2x 2 
-2xn 



Xx + r X3 



(6.5) 



and some terms that either have as a factor one of the Xi, i > 3, or are quadratic in x 2 and x 3 . 

We want to show that, near points of type 3 and some points of type 2, all focal points of order 2 are contained 
in a submanifold of codimension 3. The claim will follow if we show that the gradients of the four entries span a 
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3-dimensional space at points in U. For convenience, write r(x2, x 3 ) = olx\ + I3x 2 x 3 + 7X§. It is sufficient that the 
matrix with the partial derivatives with respect to Xi for i — 1,2,3 of the four entries have rank 3: 



1 


2 











-2 





2a 


& 


1 


P 


2 7 



The claim holds for Xi small, i > 3, unless a = and (3 = 2. This covers points of type 3. We say a point of type 2 
has type 2a if the rank of the above matrix is 3. Otherwise, the polynomial r looks: 

r(x 2 ,x 3 ) = 2x 2 £3 + 7^3 = 2x 3 (x 2 + ^x 3 ) 

We say a point of type 2 has type 2b if r has the above form and — 1 < ^ < 1. We will show that there are 
integral curves of r arbitrarily close to the one through z without focal points near z, which contradicts property 3 
in Proposition 18,31 

Take a ray t — > Cx 3 (t) passing through a point (0,0,X3,0, ... ,0). The determinant of 16.41 along the ray is: 

m - S(c(*))g(cw) - S(c(*))g(c(*)) 

= t 2 + t(Ax 2 + 27x3) + (4a| + 47x2x3 + 4x 2 ) + #3(2:3, t) 
= {t + 2x 2 + 7 x 3 ) 2 + (4 - 7 2 )x 2 + R 3 (x 3 ,t) 
>c(t 2 +x 2 )+R 3 (x 3 ,t) 

for a remainder R 3 of order 3. Thus there is a S > such that for any x 3 ^ and \t\ < 5, \x 3 \ < 5, ( X3 (t) is not a 
focal point. 

We have already dealt with points of type 3, 2a and 2b. Now we turn to the rest of points of type 2 (type 2c). 
We have either 3>lor^<— 1. We notice that x\ — x\ < iff |x2 < |x3|, but whenever |x 2 < \x 3 \, the sign 
of r(x2, x 3 ) is the sign of 7. Thus the second order part of F maps U into the complement of points with negative 
second coordinate and whose third coordinate has the opposite sign of 7. 

A similar argument as the one for type 1 points yields a contradiction with the balanced condition. If, for 
example, 7 > 2, none of the following points 

x k = (t\ -(C + l)(t k f, -{C + l){t k f, ..0, ) 

is in F(U), for t k — > 0. But then we can carry a vector other than (1,0,..., 0) as we approach F(xq). 

Type 4. Let z be a focal point of order 2. We show now that the image of the points of type 4 inside U z has 
Hausdorff dimension at most n — 3. U z is an open set around an arbitrary point z of order 2, and thus the result 
follows. 

First, we find that for any point x of type 4, we have d 2 F(v$w) = for all v,w £ ker d x F, making the 
computation in the special coordinates at x £ U z (see section [8TTI for the definition of d 2 F). 

Then we switch to the special coordinates around z. In these coordinates, the kernel of dF at x is generated by 

all and afj- Thus dx~i~ = for i>j = 2,3 at any point x £ U z of type 4. 

The set of focal points of order 2 is contained in the set H = {-q^{x) = 0}. This set is a smooth hypersurface: 

the second property in 18.31 implies that g x f* ^ at points of H. At every focal point of type 4, the kernel of dF 
is contained in the tangent to H. Thus focal points of type 4 are focal points of the restriction of F to H. The 
Morse-Sard-Federer theorem applies, and the image of the set of points of type 4 has Hausdorff dimension n — 3. 

□ 

Proof of Theorem ] 2. 11 Follows immediately from the above, setting = F(Q 2 ). □ 

7 Structure up to codimension 3 

This section contains the proof of [2. 2\ splitted into several lemmas. All of them are known for cut loci in riemannian 
manifolds, but we repeat the proof so that it applies to balanced split loci in Finsler manifolds. 

Definition 7.1. We sayp £ S is a cleave point iff R p has two elements X 1 and X 2 , with (p, X 1 ) = (F(yi), dF yi (r yi )) 
and (p,X 2 ) = (F(y2),dFy 2 (r y2 )), and both dF yi and dF y2 are non-singular. 
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Proposition 7.2. C is a (n — 1)- dimensional manifold. 



Proof. Let p = F{y\) = F(y 2 ) be a cleave point, with R p — {dF yi (r),dF y2 (r)}. We can find a small neighborhood 
U of p so that the following conditions are satisfied: 

1. U is the diffeomorphic image of neighborhoods U± and U 2 of the points y\ and y 2 . 

Thus, the two smooth vector fields X q — dF\u 1 (r) and X 2 — dF\u 2 (r) are defined in points q £ U. 

2. At all points q £ U, R q C {X^, X 2 }. Other vectors must be images of the vector r at points not in U\ or U 2 , 
and if they accumulate near p there is a subsequence converging to a vector that is neither X\ nor X 2 . 

3. Let Ti be an hypersurface in U\ passing through y\ and transversal to X±, and let T = F(T). We define local 
coordinates q — (x,t) in U, where x £T and t £ R are the unique values for which q is obtained by following 
the integral curve of X 1 that starts at x for time t. U is a cube in these coordinates. 

We will show that S is a graph in the coordinates (x, t). Let Ai be the set of points q for which R q contains X q , 
for i = 1,2. By the hypothesis, S = Aid A 2 . 

Every tangent vector v to S at q £ S (in the sense of I4.2|l . satisfies the following property (where X is the dual 
covector to a vector X £ TM.): 

X l (v) = max Y(v) 
YeR P 

which in this case amounts to X 1 ^) = X 2 (v), or 

v £ kcriX 1 - X 2 ) 

We can define in U the smooth distribution D = ker(X 1 — X 2 ). S is a closed set whose approximate tangent 
space is contained in D. 

We first claim that for all x, there is at most one time to such that (x, to) is in S. If (x, t) is in A\, R( x ,t) contains 
X 1 and, unless (x, s) is contained in A\ for s in an interval (t — e, t), we can find a sequence (x n , t n ) converging to 
(x,t) with t n y t and carrying vectors X 2 . The incoming vector is X 1 , but 

X 2 ^ 1 ) < X l {X r ) = 1 

which contradicts the balanced property. Analogously, if R( x ,t) contains X 2 there is an interval (t, t + e) such that 
(x,s) is contained in A 2 for all s in the interval. Otherwise there is a sequence (x n ,t n ) converging to (x,t) with 
t n \ t and carrying vectors X . The incoming vector is — X 1 , but 

-1 = X^-X 1 ) < ^(-X 1 ) 

which is again a contradiction. The claim follows easily. 

We show next that the set of x for which there is a t with (x,t) £ S is open and closed in T, and thus S is 
the graph of a function h over L. Take (x, t) £ U D 5 and choose a cone Z? e around Ae. We can assume the cone 
intersects dU only in the x boundary. There must be a point in S of the form (x',t') inside the cone for all x' 
sufficiently close to x: otherwise there is either a sequence (x n ,t n ) approaching [x, t) with t n > h+(x) (h being the 
upper graph of the cone D £ ) and carrying vectors X 1 or a similar sequence with t n < h-(x) and carrying vectors 
X 2 . Both options violate the balanced condition. Closedness follows trivially from the definition of S. 

Define t = h(x) whenever (a;, t) £ S. The tangent to the graph of h is given by D at every point, thus S is 
smooth and indeed an integral maximal submanifold of D. □ 

Remark. It follows from the proof above that there cannot be any balanced split locus unless D is integrable. 
This is not strange, as the sister notion of cut locus does not make sense if D is not integrable. 

We recall that the orthogonal distribution to a geodesic vector field is parallel for that vector field, so the 
distribution is integrable at one point of the geodesic if and only if it is integrable at any other point. In particular, 
if the vector field leaves a hypersurface orthogonally (which is the case for a cut locus) the distribution D (which 
is the difference of the orthogonal distributions to two geodesic vector fields) is integrable. It also follows from 12.31 
that the characteristic vector field in a Hamilton-Jacobi problem has an integrable orthogonal distribution. 
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Remark. We commented earlier on our intention of studying whether a balanced split locus is actually a cut locus. 
The proof of the above lemma showed there is a unique sheet of cleave points near a given point in a balanced split 
loci. It is not too hard to deal with the case when all incoming geodesies are non-focal, but focal geodesies pose a 
major problem. 

Proposition 7.3. The set of points p £ S where co(R*) has dimension k is (n — k)-rectifiable. 

Proof. Throughout the proof, let X be the dual covector to the vector X £ TM . 

Let p n be a sequence of points such that co (Rp n ) contains a fc-dimensional ball of radius greater than S. Suppose 
they converge to a point p and v Pn (p) converges to a vector 77. 

We take a neighborhood U of p and fix product coordinates in 7r _1 (t/) of the form U x R" . Then, we extract 
a subsequence of p n and vectors X\ £ R Pn such that X n converge to a vector X 1 in R p . Outside a ball of radius 
cS at X*, where c is a fixed constant and n >> 0, there must be vectors in Rp n , and we can extract a subsequence 
of p n and vectors X 2 converging to a vector X 2 such that X 2 is at a distance at least cS of X . Iteration of this 
process yields a converging sequence p n and k vectors 

X n , ■ ■ , X n £ Rp n 

converging to vectors 

X 1 , .., X k £ R p 

such that the distance between X k and the linear span of X 1 ,..X k ~ 1 is at least cS, so that coV* contains a 
fc-dimensional ball of radius at least c'S. 

The balanced property implies that the X 3 evaluate to the same value at 77, which is also the maximum value 
of the Z(rj) for a vector Z in R p . In other words, the convex hull of the X 3 belong to the face of R* that is exposed 
by 77. If coR* is fc-dimensional, r] belongs to 

(coi?*) = {v £ T P M : (w,v) is constant for w £ coR*} 

= (v £ T p M : (X,v) is constant for X £ 

which is a n — k dimensional subspace. 

Let S| be the set of points p £ S for which co R* is /c-dimensional and contains a fc-dimensional ball of radius 
greater than or equal to 5. We have shown that all tangent directions to S| at a point p are contained in a n — k 
dimensional subspace. We can apply theorem 3.1 in [AACj to deduce is n — k rectifiable, so their union for all 
5 > is rectifiable too. 

□ 



8 Appendix: Finsler geometry and exponential maps 

Definition 8.1. The dual one form to a vector V £ T p M with respect to a Finsler metric if is the unique one form 
lu £ T*M such that u)(V) = <p(V) 2 and ui\h = 0, where H is the hyperplane tangent to the level set 

{W £ T p M : tp(W) = tp(V)} 

atV . It coincides with the usual definition of dual one form in Riemannian geometry. 

For a vector field, the dual differential one-form is obtained by applying the above construction at every point. 



W 3 = —(P> V ) 



Remark. In coordinates, the dual one form w to the vector V is given by: 

dip 

Actually <p is 1-homogeneous, so Euler's identity yields: 

and, for a curve j(-e,e) — > T p M such that 7(0) = V, <£(7(*)) = <p(V) and 7'(0) = z, 

d 

WjZ ] = Tj£|t=o¥>(7(*)) = 
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Remark. The hypothesis on H imply that the orthogonal form to a vector is unique. 

Definition 8.2. The orthogonal hyperplane to a vector is the kernel of its dual one form. The orthogonal distribution 
to a vector field is defined pointwise. 

There are two unit vectors with a given hyperplane as orthogonal hyperplane. The first need not to be the opposite 
of the second unless H is symmetric (H(—v) — H(v) ). We can define two unit normal vectors to a hypersurface 
(the inner normal and outer normal). 



8.1 Regular exponential map 

The following proposition states some properties of a Finsler exponential map that correspond approximately to 
the definition of regular exponential map introduced in |Wj : 

Proposition 8.3. In the settina \Kl\ the following holds: 

• dF x {r x ) is a non zero vector in T F / X ^M . 

• at every point x G V there is a basis 

B = {Vl, ..,!)„} 

of T X V where r — v\ and v%,^,Vk span \erdF x , and such that: 



1 = [dF(v 1 ),d 2 F~^v 2 ), . . . d 2 F~(jfv k ),dF(v k+ i), . . . dF(v n )} 



B 

is a basis of T F ^M , where d 2 F(r$V2) is a representative ofd 2 F(r\j,V2) G T F ^M / dF{TV x ) . 

• Any point x G V has a neighborhood U such that for any ray 7 (an integral curve of r), the sum of the 
dimensions of the kernels of dF at points in 7 Pi U is constant. 

• For any two points x\ ^ X2 inV with F{x\) — F(x2), dF Xl (r Xl ) =^ dF X2 (r X2 ) 

Proof. The first three properties follow from the work of Warner |Wl Theorem 4.5] for a Finsler exponential 
map. We emphasize that they are local properties. The last one follows from the uniqueness property for second 
order ODEs. We remark that the second property implies the last one locally. Combined, they imply that the map 
p — » {F(p), dF p r) is an embedding of V into TM. 

Indeed, properties 1 and 3 are found in standard textbooks ([M])- For the convenience of the reader, we recall 
some of the notation in [W] and show the equivalence of the second property with his condition (R2) on page 577. 

• A second order tangent vector at x G R™ is a map a: 

d 2 f 



*(/) - E 



a 

1,3 



lJ dxidxj 



• The second order differential of F : V — * M at x is a map d x F :T X V T X M defined by: 

d 2 x F(a)f = a(foF) 

• The symmetric product v$w of v G T X V and w G T X V is a well defined element of T X V/T X V with a represen- 
tative given by the formula: 

(v$w)f = ~(v(w(f))+w(v(f))) 
for arbitrary extensions of v to w to vector fields near x. 

• The map d x F induces the map d 2 F : T 2 V X /TV X — ► T F ^M/dF(TV x ) by the standard procedure in linear 
algebra. 

• For x G V, v G T X V and w G ker di 7 ^, d 2 F(v$w) makes sense as a vector in the space T F ^M/dF{TV x ). For 
any extension of v and w, the vector d 2 F(v$w) is a first order vector. 

Thus, our condition is equivalent to property (R2) of Warner: 

At any point x where ker di^ 7^ 0, the map d 2 F : T 2 V X /TV X — > T 2 M F ^/dF(TV x ) sends (r^jjker dF x isomor- 
phically onto T F(x) M/dF(TV x ). 

Finally, we recall that d x F(v) is the Jacobi field of the variation F((j) t (x + sv)) at F(x), where (j) t is the flow of 
r and, whenever d x F(v) = 0, d 2 F(v$r) is represented by the derivative of the Jacobi field along the geodesic. 
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